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Abstract 

Complexity of dynamical networks can arise not only from the complexity of the 
topological structure but also from the time evolution of the topology. In this paper, we 
study the synchronous motion of coupled maps in time-varying complex networks both 
\ analytically and numerically. The temporal variation is rather general and formalized as 

being driven by a metric dynamical system. Four network models are discussed in detail in 
which the interconnections between vertices vary through time randomly. These models 
' are 1) i.i.d. sequences of random graphs with fixed wiring probability, 2) groups of graphs 

, with random switches between the individual graphs, 3) graphs with temporary random 

failures of nodes, and 4) the meet-for-dinner model where the vertices are randomly 
grouped. We show that the temporal variation and randomness of the connection topology 
can enhance synchronizability in many cases; however, there are also instances where they 
reduce synchronizability. In analytical terms, the Hajnal diameter of the coupling matrix 
sequence is presented as a measure for the synchronizability of the graph topology. In 
^ , topological terms, the decisive criterion for synchronization of coupled chaotic maps is 

QQ ' that the union of the time-varying graphs contains a spanning tree. 

\ PACS 05.45.Ra(Coupled map lattices); 05.45.Xt(Synchronization, coupled oscilla- 

^SJ ' tors); 02.50.Ey(Stochastic processes). 

(N 



1 Introduction 



Synchronization of coupled maps in networks is presently an active research topic [1]. It 
represents a mathematical framework that on the one hand can elucidate - desired or 
. undesired - synchronization phenomena in diverse applications. On the other hand, the 

^ ' synchronization paradigm is formulated in such a manner that powerful mathematical 

techniques from dynamical systems and graph theory can be utilized. A standard version 
of the network of coupled maps, coming from the well-known coupled map lattices (CML) 
[2], can be formalized as follows: 



x^{t + l)^ f{x\t))+y^U,f{x\t)),i^l,--- (1) 
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where t G Z+ = {0, 1,2,-- - ,}, x^{t) is the state variable of vertex j, / : R — *■ R is 
a differentiable map, and L = [Lij]"''^^^ G jjmxm -g diffusion matrix, which is 
determined by the topological structure of the network and satisfies Lij > for all 
i ^ j and E^i -^y = for all i = I,-- - , m. Let x = x^, . . . , x™]^ e R™, 
F(a;) = [/(a;^), /(a;^), . . . , /(x™)]^, and G — I„i + L, where denotes the identity 
matrix of dimension m. Then, Eq. ([T]) can be rewritten in the matrix form: 

x{t + I) = GF{x(t)) (2) 

where G — [Gij]™ ^;^ G jjmxm (^gj^otes the coupHng and satisfies > for i j and 
SjLi Gjj = 1 for alH = 1, • • • , m. Thus, if Gu > holds for alH = 1, • • • , m, then G is 
a stochastic matrix. 

This dynamical system formulation contains two aspects. One of them is the reaction 
dynamics at each vertex of the network. The other is the coupling structure, that is, 
whether and how strongly, the dynamics at one vertex is directly influenced by the states 
of the other vertices. This influence can be described by notions of graph theory. Hence, 
the coupling matrix G corresponds to a graph T = [V, £], where V = {1, 2, • ■ • , m} denotes 
the vertex set and £ = {e^ } denotes the edge set such that there exists a directed edge 
from vertex j to vertex i if and only if Gy > 0. 

Synchronous dynamics in complex networks have recently attracted increasing atten- 
tion [1,3-7]. Linear stability analysis was used and transverse Lyapunov exponents were 
introduced to analyze the influence of the topological structure of networks [6]. Ref. [7] 
has related the ability to synchronize chaotic maps to the existence of a spanning tree 
in the corresponding graph. However, synchronization analysis has so far been mostly 
limited to autonomous systems, where the interactions between the state components are 
static. In [6], a generalized criterion guaranteeing synchronization in the model ([2]) is 
proposed as follows: 

log|Ai|+^<0, (3) 

where is the Lyapunov exponent of the uncoupled system s{t + 1) = f{s{t)) and Ai is 
the eigenvalue of the coupling matrix G with the second largest modulus, noting that the 
largest eigenvalue has a modulus of 1. 

Many real-world applications from the social, natural, and engineering disciplines 
include a temporal variation of topology of the network. In communication networks, 
for example, one must consider dynamical networks of moving agents. In this case, 
some of the existing connections can fail simply due to occurrence of an obstacle between 
agents [8]. Also, some new connections may be created when one agent enters the effective 
region of other agents [9-11]. Furthermore, this temporal variation of topology involves 
randomness. In [8-10], consensus in multi-agent networks was considered where the state 
of each vertex is updated according to the states of its connected neighbors with switching 
connecting topologies. The consensus protocol of multi-agent dynamical networks can 
generally be formalized in discrete-time form as 

m 

x*(i + l) =^Gy(t)x^(t), z = l,---,m, (4) 

where [Gij{t)]™j^^, t G Z+, are stochastic matrices. It was proved in Ref. [11] that 
the connectivity of the switching graphs plays a key role in the consensus dynamics of 
multi-agent networks with switching topologies. Some papers from the recent literature 
[12] studied synchronization of continuous-time dynamical networks with time-varying 
topologies; however, the time- varying couplings were specific, with either symmetry, node 
balance, or fixed time average. 



INTRODUCTION 



In this paper, we study the local complete synchronization of networks of coupled 
maps with time- varying couplings: 

xit + l) = G{e'''^uj)F{x{t)). (5) 

Here, 6*^*^- represents a metric dynamical system {Q,T,P,9^^^}, where fl is the state 
space, is the cr-algebra, P is the probability measure, and 6'*^*^ is the semi-flow satis- 
fying 6l(*+'') = 6l(*) o 6'("), where 9^°'> is the identity map, G(6'(*)w) = [G»j(6'(*)u;)]™-^i € 
]^mxm (denotes the coupling matrix at time t and is measurable on and F{x) — 

[/(xi), • • • , f{xn)]^ is a differentiable function. 

Thus, Eq. ([5|) is a random dynamical system. For more details on random dynamical 
systems, we refer to the textbooks [16]. This form of time- varying coupling is rather 
general and includes the deterministic case, where G(-) can be regarded as a known 
function of time t, as well as the stochastic case, where G(-) can be regarded as being 
enforced by a stochastic process {C*}tez+: namely, G(^*). 

Accordingly, we denote time varying graphs by {r(0(*)tj)}fgz+ ■ Define r(0(*)w) = 
[V,£{9'^*^iL!)], where V = {1,2,- •■ ,m} denotes the fixed vertex set and £"(0^*^^^) — 
{cij {d'-*^ Lj)} denotes the edge set of the graph at time t, i.e., edge eij{d'-*^uj) exists if 
and only if Gij{9'^*''^uj) > 0. So, the coupling matrix GiO'^^'^uj) might be a function of the 
coupling graph topology. 

Local complete synchronization (synchronization for short) is defined in the sense that 
the differences between states of vertices of the coupled dynamical system ([S|) converge 
to zero whenever the initial state of each vertex is picked sufficiently near the attractor 
of the uncoupled system and their differences are sufficiently small, i.e., 

lim \\x\t) -x^{t)\\ = 0, i,i = I,- - • ,m. (6) 

t — 'OC 

For a more geometric definition, suppose that the uncoupled system s{t -t- 1) = f{s{t)) 
possesses an attractor (see Ref. [17] for details), which we denote by A. Define 

S = {[x\x^,--- ,x"']'^ £ M™ : x' ^x^, = 1,.. . ,m} 

which is an invariant subspace of Eq. ^ . Let A™ denote the Cartesian product Ax - ■ - xA 
(to times). We define the synchronization manifold by the set A = SO A™ — {[x, - ■ ■ ,x] : 
X £ A}. In this sense, synchronization is equivalent to the stability of A. 

The purpose of this paper is to study the synchronization of the coupled map network 
Q with time-varying topology. Here, the topology is generally supposed to be driven by 
a metric dynamical system and the coupled network can be regarded as a random dynam- 
ical system. We present sufficient conditions guaranteeing synchronization. Furthermore, 
we show that the property that the union of the time- varying graphs contains a spanning 
tree is very important for the network's ability to synchronize chaotic maps. Addition- 
ally, we present several time-varying network models and study the synchronization of 
coupled maps on these dynamical networks. The topological structures of these models 
vary in time and include randomness. Generally, the collections of interconnections in 
these networks can be regarded as Markov chains. Besides illustrating the theoretical 
results, we also focus on the variation of synchronizability of each model, which is quanti- 
tatively measured with respect to several parameters in the model. As we show, temporal 
variation and randomness can enhance synchronization in some cases. Further examples 
indicate that the communication between vertices in the dynamical networks might play 
an important role in synchronizability. 



THEORETICAL ANALYSIS 



2 Theoretical analysis 

In this section, we present theoretical results on synchronization of coupled map networks 
with time-varying couplings. The mathematical results have been proven in detail in our 
companion papers [14, 15]. Our main tool to investigate the synchronous motion of the 
coupled system ([S]) is the Hajnal diameter, which was first introduced in Ref. [13] to 
describe the compression rate of a stochastic matrix and is defined as follows: 



Definition 1 For a matrix G with row vectors gi , ■ 
the Hajnal diameter of G is defined as 

diam(G') — max \\qi — qj \\. 



, Qm and a vector norm || • || in . 



(7) 



From this definition, synchronization of the coupled system ^ can equivalently be stated 



as 



lim diam([a;i(t), • • • , Xm{ty\^) = 0. 



(8) 



We can extend this concept to matrix sequences driven by a dynamical system: Q{lij) = 
{G{9^^^u)}t>o : ^ 2* " for any w g 17, where 2* ' denotes the set composed of all 
subsets of M™'™. For a matrix sequence Q, its Hajnal diameter at initial data w G is 
defined by 



r t-i 



diam((y(a;)) = lim {diam 



k=0 



(9) 



where Y\ denotes the left matrix product: nfc=i ~ An x An-i x • • • x Ai. One can 
see that diam(Cy(a;)) < implies that the differences between rows of the infinite matrix 
product Ht^o G{9^*''uj) converge to zero as t goes to infinity. 

Let s{t) be the synchronized state solution satisfying s{t + 1) = f{s{t)) for all t > 0. 
Let Sx{t) = x{t) — s{t). Linearizing the system (O about s{t) gives 



6x{t + 1) = f{s{t))G{e'-*^uj)Sx{t). 



(10) 



Note that 



diam 



- k=0 
r t-1 

diam 



nG(0^'^c.)/'(/W(so)) 



fc=0 



1=0 



(11) 



Then, the Hajnal diameter of the variational system (fTU]) equals to diam(t/(aj))e^, where 
fj, denotes the maximum Lyapunov exponent of the attractor A of the uncoupled system. 



1 *-i 

maxlim-^log|/(/W(so))|. 



k=0 



This leads to the following condition 

diam(e;(a)))e^ < 1, 



(12) 



(13) 



which guarantees that the variable vector x{t) can be synchronized by picking the initial 
data of 6''*'- as ui. 
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Similar to the case of static network topology, we can extend the transverse Lyapunov 
exponent for the matrix sequence Q in direction v £ M™ as: 



a{Q,uj,v) = lim - log 

t — ^oo t 



1 



(14) 



Along the synchronization direction eo = [1, 1, • • • , 1]^, 

one has fj{Q ,uj,eQ) — since G{-) has a common row sum of unity. Let = cto > cri > 
> • • • > cr„i be the Lyapunov exponents for the initial condition w, counted with 
multiplicities. We have ai{u)) = logdiam(C/(a;)) according to lemma 2.7 in Ref. [14]. 
Then, the condition (1131) can be rewritten as 



cri+/i<0. (15) 
If is satisfied, then the coupled system ([5]) can synchronize. 

Remark 1 By Proposition 4-4 R^f- [14]> one can see that the criterion (0) for static 
networks is a direct consequence of the criterion hl5\). 

We apply the above results to the case where the time- varying coupling is induced by a 
homogeneous Markov chain {(T*}tgz+ defined on a finite state space with an irreducible 
transition probability matrix. Also, a homogeneous Markov chain can be regarded as 
a dynamical system (fl, JF, Ptt, 6''*''-) as described in the appendix. We now consider a 
coupled map network with Markov jump topologies: 

m 

x\t + l)=Y,G,,{a')f{x^{t)), z = l,...,m (16) 



or in matrix form: 



x{t + l) = G{a')F{x{t)). (17) 



Results in Ref. [15] indicate that in this case, logdiam(^(a;)) = cri(w) exists and is a non- 
random number for almost every w £ $7. Hence for simplicity we can write diam(f/(cLi)) 
as diam(C?) and cri(w) as cxi. From (|13p. one can obtain the criterion for synchronization 
of coupled maps (fT7|) as 

logdiam(g) +^ < 0. (18) 

According to the equivalence, we can rewrite the condition (jl8p as the inequality p5p . 
From the criteria (fT5|) -(fT5 |) . the Hajnal diameter diam(C/), or equivalently, ai, can be used 
to measure the synchronizability of a Markov jump graph topology process. The ques- 
tion then arises under what conditions this graph process can synchronize some chaotic 
dynamics, i.e., when does it hold that diam(C7) < 1- The following result comes from 
Theorem 4.2 in Ref. [14] and the theory of Markov chains [21]. 

Theorem 1 Suppose that G(-) has all diagonal elements positive and the transition prob- 
ability matrixT is irreducible. Then, diam(Q) < 1 if and only if the graph union IJjgjv ^i^) 
possesses a spanning tree. 

For the detailed proof, we refer to Ref. [15]. This theorem shows that there exist cases 
when a Markov jump graph process can synchronize a chaotic map (with /i > 0) even 
though at each instant the network may be disconnected, as long as the union graph has 
a spanning tree. 
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3 Applications 

In the following, we will study the synchronous dynamics in four time-varying graph 
process models. In each model, the number of vertices is constant in time but the inter- 
connections between vertices vary, and the variation of interactions can be regarded as 
a Markov chain. We expect, on the one hand, to illustrate the theoretical results of the 
previous section, and on the other hand, to numerically analyze the synchronizability as 
computed by the largest nonzero Lyapunov exponent ai, by observing the variations of 
CTi with respect to several parameters in the models. 

The map / is chosen here as the logistic map: / = ax{l — x). We take the parameter 
a = 3.90 throughout this section (hence with the Lyapunov exponent fi w 0.5). Thus, 
we can focus on the influence of the time-varying coupling on synchronous motions by 
fixing the parameter of the coupled map, which fixes the Lyapunov exponent ^ of the 
uncoupled system. (Note that the theoretical results presented above do not depend on 
this particular choice of chaotic dynamics.) 

We realize the coupled networks via two types of coupling configurations. The first 
system is the coupled map lattice via a time varying graph process: 

m 

/(a;*(t)) + ^^A,,(i)[/(x^W-/(^Xt)], iffc.(t)>0, ^^g^ 
fix'it)), ' iffc,(t) = 0, 

where i = 1, • ■ • , m, e > is the coupling strength, A{t) denotes the adjacency matrix of 
the graph at time t, and ki{t) — (*) denotes the (in-)degree of vertex i at time 

t. Synchronization is measured by the time average of the variance of the states over the 
network: 

, rn 

K=( J2[:,^(^t)-x{t)]'' 

\ TO — 1 — ' 

^ i=l 

where x — (1/™) X]"=i ^^'^ (') denotes the time average. One can regard K as a 
function of the couphng strength e. Let [G{t)]ij — Sij{l — e) + ek^^{t)[A{t)]ij if ki{t) > 
and [G{t)]ij — 6ij otherwise, where Sij are the elements of the identity matrix /,„. Then, 
the second largest Lyapunov exponent cri of the stochastic matrix series {G(t)}tg2+ is 
also a function of e. We also define W = ai + fi, which is the largest Lyapunov exponent 
of the system (jl9p in directions transverse to the synchronization manifold. 

The second system is the dynamical multi-agent system with the logistic output func- 
tion / given above. At each vertex i, the state is the average of the values /(x*(i)) of all 
its neighbors and itself, i.e.. 



x\t + l) 



(20) 



where Afi{t) denotes the neighborhood of vertex i in graph T{t) and ki{t) is the degree of 
the vertex i at time t. So, Gij{t) = \/{ki{t) + 1) in the form ^ if vertex j is linked to 
vertex i at time t\ otherwise, Gij = 0. According to the criterion psp . the quantity ai can 
be utilized to measure synchronizability of the time varying graph process of the coupled 
system pO|) . A smaller value of ai indicates better synchronizability. The simulation 
time length is 1000 in all cases. 

3.1 I.i.d. random graphs 

In the independent-identical-distribution (i.i.d.) random graph, the edge for each pair of 
vertices can disappear or appear randomly, independent of time and other pairs of vertices 
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and following an identical distribution. This is a special case of the model introduced in 
Ref. [9]. As a realization in the present paper, at each time t, T{t) is a p-random graph 
following the famous Erdos-Renyi model [22]: for every pair we randomly put an 

edge between them with probability p and the selection is statistically independent for 
different times t and other pairs of vertices. 

We realize the coupled map networks ([T^ and (|20p in this model. Figure [T] (a) in- 
dicates that the parameter range of the coupling strength e for which synchronization 
occurs coincides with the range where W{e) < 0. This verifies the criterion (fTSll . From 
the criterion the synchronizability measure for static networks is log |Ai|, which has 
been studied in e.g. Ref. [27]. From figure [2] (a) , we observe the variation of ci with re- 
spect to p and compare it to the logarithm of the second largest eigenvalue (in modulus) 
of the coupling matrix of the static random graph in the coupled model One can see 
that the synchronizability of i.i.d. random graphs increases with increasing probability 
parameter p and is clearly better than a static random graph of the same size and with 
the same wiring probability p. This implies that in a random network, temporal varia- 
tion and randomness can increase synchronizability. Furthermore, as one would expect, 
synchronizability increases with the wiring probability p. 

3.2 Randomly switching topologies 

Randomly switching topologies were introduced in Ref. [8]. That is, the graph topology 
at time t is randomly picked from a given finite set of topologies that follows an identical 
time-independent distribution. Here, we consider two pairs of graphs, see {Ti, and 
{Fa, in figure[31 The random switch occurs between the two graphs of each pair. The 
switching signal is driven by a Bernoulli random variable v. For some constant p G (0, 1), 
if V < p then the first graph in each pair is chosen as the coupling topology; otherwise, 
the second graph is chosen. 

From figure [T] (b) , one can see that the parameter region for which ai + fi < equals 
to the region where K 0, which verifies the criterion From figure [5] (b), one can 

see that for the graph pair {ri,r2}, the synchronizability of random switching measured 
by ai is worse than either of the individual graphs (noting that the synchronizability of 
each graph can be found at the endpoints p = and p = 1). In contrast, for the graph 
pair {Fa, F4}, the synchronizability obtained by random switching is better than those of 
the individual graphs. That is to say, there exist instances where temporal variation of 
the network topology can increase or decrease synchronizability. 

3.3 Random errors 

In this model, we consider a network with random errors. If an error occurs at a vertex, 
then all connections of this vertex disappear. This model is characterized by two kinds of 
errors. One is called failure, which happens to vertices following the uniform distribution; 
the other is called attack, which happens to vertices following a selective distribution 
according to a certain statistical property of the vertices. As shown in Ref. [23], for a class 
of complex networks with inhomogeneous degree distribution (for example, the Barabasi- 
Albert (BA) model), the statistics such as shortest-path diameter and clustering can 
have good error tolerance if the errors occur as failure but they are extremely vulnerable 
for attacks based on highest degrees. As shown in Ref. [25], the synchronizability of a 
network measured by the eigenratio of the corresponding Laplacian almost does not vary 
if a vertex is randomly removed but dramatically varies by the selective removal of one 
vertex. In the present paper, we realize attack according to the connection degree of each 
vertex. Namely, errors happen merely on vertices with highest degrees. In addition, we 
add a recovery phase: Every malfunctioned vertex will recover, i.e., all its connections 
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Figure 1: Variation of K and W with respect to e. Figures (a)-(d) are plotted for coupled 
logistic maps (|19p . (a) i.i.d. random network with 200 vertices and p = 0.1; (b) networks 
with switching topologies between and r4 with switching probability p = 0.5; (c) random 
error model: beginning with a scale-free network introduced in [24] with 200 vertices and 
average degree 20, failure occurs with probability p = 0.01 and the recovery time T = 3; (d) 
meet-for-dinner model with = 200 members and subgroups of size n = 5 . In all cases, K 
is shown by solid lines (— ) and W is shown by dotted lines (— ■ — ). 
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Figure 2: Synchronizability of different graph processes. Subfigures (a)-(d) are plotted for 
the variation of synchronizabihty ai with respect to the parameters of the coupled network 
()20p . (a) The variation of ai of an i.i.d. random network with respect to the parameter p and 
log|Ai| where Ai is the second largest eigenvalue of the coupling matrix of a static random 
graph with the same p in the model ()20p . The network size is 1024. (b) The variation of cJi of 
a randomly switching network with respect to p, for the first group {ri,r2} and the second 
group {r3,r4} of Figure O (c) The variation of cxi of malfunction-and-recovery networks 
with respect to malfunction fraction p, with recovery time T = 5, for failure and attack on 
scale-free and random networks. The initial scale-free network has size N = 1024 and average 
degree 20, and the random network has N = 1024 and average degree 152. (d) The variation 
of cTi of the meet-for-dinner model with respect to n, in a network of size N = 1024. 
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will appear, after a fixed time period. We denote by p the fraction of error vertices in 
the whole vertex set, i.e., there are [N x p\ error vertices, where [-J denotes the floor 
function and N is the size of the whole network. 

Comparing the regions of the coupling strength where K ^ and < in figure [1] 
(c), one can similarly see that the inequality (jlSp can precisely predict synchronization. 
Figure[2](c) indicates the variation of the synchronizability with error occurrence. We use 
two network models, namely the Barabasi- Albert (BA) network introduced in Ref. [24] 
as a scale-free network (which has a power- law degree distribution P{k) ^ with 

7 = 3 independent of the size of the network in case of sufficiently large network sizes), 
and a random network with exponential tails introduced in Ref. [22]. One can see that 
for a random network with high degrees, owing to the homogeneity of the network, there 
is no substantial difference in synchronizability whether the malfunctioned vertices are 
selected randomly or in decreasing order of connection degree. On the other hand, a 
drastically different behavior is observed for the scale-free network. If the vertices with 
higher degrees are attacked, the synchronizability is much reduced compared to the case 
with random failure. Due to the degree inhomogeneity of the BA networks, the vertices 
with a high degree play a more important role in synchronization than those with smaller 
degrees. 

3.4 Meet for dinner 

In the meet-for-dinner model introduced in [26], a group of friends decide to meet for 
dinner at a particular restaurant but fail to specify a precise time. On the afternoon of 
the dinner appointment, they need to find a solution to decide on the meeting time. A 
centralized solution is to have an advanced conference for the whole group; however, if 
this option is unavailable, then a decentralized solution is that one meets, one at a time, 
a subset in the subgroup to collect the information of this subgroup about their expected 
meeting time, and update with this information until obtaining consensus. Here, we set 
up the model as follows. The whole group has N members. At each time interval, the 
group is randomly divided into subgroups with n members (if iV 7^ mod n, then we 
put the remaining ones into the last subgroup) and each subgraph is a complete graph. 
Furthermore, every division is stochastically independent of each other. 

The region of the coupling strength for synchronization coincides with the range where 
< in figure [1] (d). (The tiny region of apparent discrepancy near e ~ 0.45 is an 
artifact of plotting the curves with finite data points.) Interestingly, the meet-for-dinner 
model can synchronize a chaotic map / despite the fact that the graph is disconnected at 
any time. For a static disconnected graph, there exist several vertices whose dynamical 
information never reaches the others; so, obviously, a chaotic map cannot be synchronized 
by a disconnected graph. However, if the graph topology is time-varying, despite the 
disconnectedness of the network at each time, the dynamical information can reach others 
in a certain time period. Therefore, in this sense, theorem [T] implies that in some cases, 
temporal variation of the network topology can enhance synchronization. Figure [Hd) 
shows that the synchronizability of the meet-for-dinner model increases with size n of the 
subgroups. 

4 Conclusions 

In conclusion, we have presented an effective method based on the extended Hajnal di- 
ameter for matrix sequences to study the synchronization in networks of coupled maps 
with time varying topologies. As shown by the sufficient criteria guaranteeing synchro- 
nization, the Hajnal diameter of the coupling matrix sequence can be utilized to measure 
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Figure 3: Two groups of graphs {ri,r2} and {r3,r4}. Switch occurs in either group ran- 
domly. Namely, at each time, with probability p the graph topology is selected as Fi (respec- 
tively, Fs) and with probability 1 — p as F2 (resp., F4). 



network synchronizability. As shown in Sec. 3, synchronizability varies with respect to 
several parameters in time-varying network models. An intuitive interpretation is that 
the time-cost of communication between vertices might play a key role for synchronization 
of a dynamical network. The vertices in the i.i.d. random graphs have a higher chance to 
access others than in a static random graph. Attack to a network with a power-law degree 
distribution is more likely to interrupt the communication between vertices than random 
failures. However, for a random network with high average degree, attack and failure can 
cause almost equal damage in communication between vertices. When the network size 
increases, the indirect communication of two vertices can be enhanced by the time- varying 
connection structure, which can increase synchronizability. These phenomena imply that 
in some cases time-variance and randomness can enhance synchronizability. However, as 
shown in figure [2] (b), it is also possible to have decreased synchronizability. This issue 
deserves further investigation in the future. 



Appendix: Homogeneous Markov chain with finite state 
space 

A homogeneous Markov chain with finite state space and an irreducible probability tran- 
sition matrix can be regarded as a metric dynamical system with invariant probability 
{fl, J-, P,6'^*^ ■}. Its state space fl = is composed of all sequences: u = {o'*}t>o; 

its Borel cr-algebra J- = i3^^ , where B denotes all subsets of N_, has a basis of the form 
{cT*i e Bi, • • • , cr*-- G Br} for some ti < t2 < ■ ■ ■ < U and Bi e B for ah / = 1, • ■ • , r; 6* 
denotes the shift map, 9u} = {cr^*^}t>i; P^- denotes the probabihty measure induced by 
the unique invariant class of the transition probability matrix tt, which is given by 

N 

P^{a*^ eBi,--- ,a^- eBr)^ ^ ^ Tr^tni^ti^zs ■ ■ ■ tit,._,u, (21) 
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where T = {h, - ■ ■ ,tr}, and invariant through 9^*^-. Induced by different initial distri- 
butions ^, this system can have difFcrcnt probability measures P^, but they all are not 
invariant over If the invariant probability tt is ergodic in the sense that each tt/c > 0, 
then for any initial distribution ^, is absolutely continuous with P^, i.e, P| <C Ptt, which 
implies that any characteristic in the P^r almost sure sense certainly holds in the almost 
sure sense for any Pj if tt is ergodic, or equivalently, if the transition probability matrix T 
is irreducible. In this paper, we only focus on the probability measure Ptt and simplify P^- 
almost surely by "almost surely" unless denoted otherwise. By the multiplicative ergodic 
theorem for random dynamical systems [16], the multiplicative Lyapunov exponents for 
the infinite matrix sequence Hf^o ^i'^*) exist and are non-random almost surely. 
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